#include "PowerSpectrum.h"

std::ostream* dbgout=0;
bool XDEBUG = false;

int main()
{
  Universe u(0.3,0.7,0.04,0.7,0.7,0.96,-1.0,0.0);
  double z1 = 0.545;
  double z2 = 0.716;
  std::cout<<"w("<<z1<<") = "<<u.w_z(z1)<<std::endl;
  std::cout<<"w("<<z2<<") = "<<u.w_z(z2)<<std::endl;
  return 0;
}
